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METAL-POOR GLOBULAR CLUSTERS AND THE FORMATION OF THEIR HOST GALAXIES 
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ABSTRACT 

We have determined the total numbers and specific frequencies of blue, metal-poor globular clusters 
(GCs) in eight spiral and early-type galaxies. These data, along with five measurements from the 
literature, show a trend of increasing blue GC specific frequency with increasing mass of the host 
galaxy. The increase is not accounted for in a simple galaxy formation model in which ellipticals 
and their GC systems are formed by the merger of typical spiral galaxies. The data appear broadly 
consistent with hierarchical formation scenarios in which metal-poor GCs are formed over a finite 
period in the early Universe during the initial stages of galaxy assembly. In this picture, the observed 
trend is related to biasing, in the sense that the more massive galaxies of today began assembling earlier 
and therefore formed relatively more GCs during this early epoch of metal-poor GC formation. We 
discuss how comparisons of the observed specific frequency of metal-poor GCs with model calculations 
can constrain the formation redshift of these objects. 

Subject headings: galaxies: star clusters — galaxies: formation — galaxies: spiral — galaxies: elliptical 
and lenticular, cD 



1. INTRODUCTION 

The precise way that galaxies like our own Milky Way 
form is one of the central issues in astrophysics and cos- 
mology today. One method for investigating a galaxy's 
origins is to study its stellar populations and use their 
properties to piece together a coherent picture of the dif- 
ferent stages of the galaxy's history. Globular clusters 
(GCs) make particularly useful probes of galaxy forma- 
tion because these luminous, compact star clusters are 
easy to see at large distances and appear in all types of 
galaxies, and the oldest GCs have ages approaching the 
age of the Universe (e.g.. Ashman & Zepf 1998, Chaboyer 
et al. 1998). Furthermore, GCs are thought to form dur- 
ing periods of intense star formation and therefore trace 
the major star-forming episodes in a galaxy's past. Care- 
fully quantifying the observable properties of the entire 
system of GCs in a galaxy leads to important insights 
into that galaxy's origins. 

With this rationale in mind, we initiated a ground- 
based imaging survey of the GC systems of several gi- 
ant galaxies (Rhode & Zepf 2001, 2003, 2004; hereafter, 
RZOl, RZ03, RZ04). The survey was designed to produce 
reliable determinations of the total number of GCs in 
each galaxy and to quantify the spatial and color distri- 
butions of the GCs out to large galactocentric radius. We 
use large-format CCDs to observe the entire radial extent 
of the GC systems. Furthermore, our survey methods 
(e.g., three-color photometry) facilitate the separation of 
real GCs from contaminating foreground and background 
sources. As a result, we have significantly reduced the 
levels of contamination in our samples and reduced the 
errors on our calculated GC specific frequencies by a fac- 
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tor of 2-3 compared to past studies (RZOl, RZ03, RZ04). 

The survey's main purpose is to help answer questions 
about how galaxies and their GC systems form. In par- 
ticular, we aimed to test the predictions of models for 
the origins of ellipticals. Two important properties of 
elliptical galaxy GC systems are that they are typically 
more populous than those of spiral galaxies and that they 
often show two peaks in their color or metallicity distri- 
butions (Ashman & Zepf 1998 and references therein). 
Ashman & Zepf (1992; hereafter AZ92) predicted this 
bimodality when they proposed that ellipticals and their 
GC systems are created in mergers of disk galaxies. In 
this picture, a metal-poor (blue) GC population in the 
elliptical comes from the progenitor spiral galaxies, and a 
metal-rich (red) population is formed during the merger. 
A consequence of this scenario is that the relative num- 
bers of metal-poor GCs should be similar in spiral and 
elliptical galaxies. This can only be tested with a careful 
accounting of the total numbers and specific frequencies 
of GCs (both blue and red) in a sample of early and late- 
type giant galaxies. Such data are also useful for testing 
other models proposed to account for the observed bi- 
modal metallicity distributions (e.g., Forbes, Brodie & 
Grillmair 1997, Cote, Marzke, & West 1998). 

More broadly, the low metallicities (the mean [Fe/H] 
for the metal-poor Milky Way GCs is —1.6; Armandroff 
& Zinn 1988) and extended spatial distributions (RZOl, 
RZ04) of metal-poor GCs suggest that they were some 
of the first structures to form in the Universe. There- 
fore we can use what we learn about metal-poor GC sys- 
tems to probe structure formation at early epochs. In 
hierarchical cosmological models, galaxy halos that are 
more massive and are located within dense environments 
will tend to have a higher proportion of their mass ac- 
cumulated by a given redshift compared to less-massive 
galaxies in poorer environments (e.g., Sheth 2003). The 
assumption of a preferred formation epoch and constant 
formation efficiency for metal-poor GCs then leads to the 
prediction that higher-mass galaxies will have larger rela- 
tive numbers of metal-poor GCs. Furthermore, the exact 
spatial distribution of the metal-poor GCs out to large 
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radius constrains the formation redshift and the masses 
of the structures within which they formed (Santos 2003 
and references therein). 

We have recently completed the analysis for three more 
galaxies from the survey (Rhode & Zepf 2005) bring- 
ing our total sample to eight spiral, elliptical, and SO 
galaxies. We can combine these results with those from 
other studies to address some of the questions that mo- 
tivated the survey and to investigate the possible cosmo- 
logical bias of metal-poor GCs. In this letter, we present 
our results for the specific frequency of metal-poor GCs 
and discuss what they might imply about how and when 
galaxies and their GC systems formed. 

2. SPECIFIC FREQUENCY OF METAL-POOR GLOBULAR 
CLUSTERS 

2.1. Survey Data and Analysis 

Details of our data analysis methods are given in a 
series of previous papers (RZOl, RZ03, RZ04). Briefly, 
to estimate the number of metal-poor GCs in a target 
galaxy, we begin by imaging the galaxy in three broad- 
band optical filters {BVR or BVI). GC candidates are 
selected by their radial profiles (GCs appear unresolved 
from the ground at the 10—20 Mpc distances of our 
targets), apparent magnitudes, and colors. The spatial 
distribution of GC candidates for each galaxy is then 
constructed and corrected for incompleteness due to the 
magnitude depth of the images, missing spatial coverage 
around the galaxy, and contamination from non-GCs. 
Integrating this distribution out to the radius at which 
the GC surface density drops to zero yields an estimate 
of the total number of GCs in the system. 

As mentioned earlier, the broadband color distribu- 
tions of many giant galaxies often show two peaks, indi- 
cating the presence of two populations of GCs (e.g., Zepf 
& Ashman 1993; Gebhardt & Kissler-Patig 1999; Kundu 
& Whitmore 2001; RZOl; RZ04). The color difference 
represents a difference in the metallicity of the two pop- 
ulations in nearly all cases, because the GC systems are 
typically more than a few Gyr old (Zepf & Ashman 1993; 
Ashman & Zepf 1998). In a few well-studied galaxies, the 
spatial distributions and kinematics of the blue (metal- 
poor) and red (metal-rich) GCs have been shown to differ 
from each other, supporting the interpretation that these 
are distinct populations created in different star forma- 
tion episodes (e.g., Zinn 1985, Zepf et al. 2000). To cal- 
culate a total number of blue GCs in our galaxies, we 
use the KMM mixture-modeling code (Ashman, Bird, & 
Zepf 1994), which (in addition to testing whether a sin- 
gle or multiple Gaussian function provide a better fit to 
a distribution) locates the peaks and dispersions of the 
Gaussians that best fit a distribution and provides an es- 
timate of the fractional contribution of each population. 
The code is run on a sample of GC candidates for which 
our detection is complete in all three filters. 

For the early-type galaxies in the sample, the sepa- 
ration between the blue and red GCs occurs aX B — 
R ~1.2-1.3, and 60-70% of the GCs are blue (RZ04). 
The spiral galaxies in the study often have poorer statis- 
tics, making definitive separation of the candidates into 
blue and red populations difficult or impossible. In such 
cases we counted as metal-poor the GCs bluer than 
B — R^ 1.23, the typical location of the separation in the 
color distributions of the early- type galaxies (RZ04). The 



gap between the metal-poor and metal-ri ch GCs in the 
Milky Way occurs at B-R - 1.25 l|Harris 1 1996). indicat- 
ing this is a reasonable choice. The estimated blue/red 
proportions of GCs in the spiral galaxies were not sensi- 
tive to the exact location of the B ~ R separation; using 
1.23 or 1.25 produced the same results within the errors. 
The total number of blue GCs was computed for each 
galaxy by multiplying the estimated blue fraction by the 
total number in the system. 

2.2. Measurements from the Literature 

To supplement our data, we searched the literature for 
well-determined measurements of the numbers of blue 
GCs in other giant galaxies. We used the following cri- 
teria for inclusion of a particular measurement here: at 
least 50% of the estimated radial extent of the GC system 
must have been observed; data must have been acquired 
in at least two filters, so the GC color distribution could 
to be quantified; and estimates of the total number of 
GCs and the blue fraction were given, or could be de- 
rived in a fairly straightforward way from the published 
data. Finally, the measurement was included only if its 
l-cr error was ^40%. A comprehensive search of refer- 
eed journal papers yielded measurements that met these 
criteria for three elliptical galaxies. We also include esti- 
mates of the number of blue GCs in the two most well- 
studied giant galaxies, the Milky Way and M31, in our 
analysis. The data for all the galaxies and the references 
from which they are derived are given in Tabled 

2.3. Results 

In order to compare the blue GC populations of the 
galaxies, we have calculated the mass-normalized number 
of blue GCs, or T biuc, for each. This is a variation of the 
T parameter introduced bv iZeof fc Ashman . (,1993.'! and 
is defined as: 

_ iVGc(blue) 

- Mg/109 Mq ^ ' 

where Nqc (blue) is the number of blue GCs and Mq is 
the stellar mass of the host galaxy. Normalizing Nqc 
by galaxy mass facilitates the comparison of the GC sys- 
tems of galaxies with a range of morphological types. 
Normalizing instead by galaxy luminosity within a given 
waveband — e.g., as in the s pecific frequency ^jv of 
iHarris fc van den Bergh I l)198lD — makes it difficult to 
directly compare numbers of GCs in spirals and ellipti- 
cals with very different stellar populations. 

Figure n shows Tbiuc for 13 giant galaxies: eight from 
our survey and five from other work. To calculate galaxy 
masses, we assume mass-to-lig ht ratios {M/Ly) follow- 
ing those of lZepf fc Ashman I 1^93): M/Ly = 10 for 
ellipticals, M/Ly = 7.6 for SO galaxies, and M/Ly be- 
tween 5.0 and 6.1 for spirals, depending on Hubble type. 
The relevant quantity for comparing galaxies across Hub- 
ble types is the relative mass-to-light ratio as a function 
of type, and analyses of both stellar populations and dy- 
namics yield relative mass-to-light ratios similar to those 
adopted here (e.g.. Ashman & Zepf 1998). Kauffmann 
et al. (2003) analyzed Sloan Digital Sky Survey data and 
also found similar or slightly larger mass-to-light ratio 
differences among galaxy types. 
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Fig. 1. — Mass-normalized numbers of blue GCs for 13 giant 
galaxies versus the log of the galaxy stellar mass. Circles are cluster 
ellipticals, squares are early- type field galaxies, and triangles are 
field spiral galaxies. Filled symbols are used for our data and the 
Local Group galaxies; open symbols are data from other studies. 



Tbiuc is shown in the figure with difi^erent symbols to 
indicate galaxy type and environment: circles for clus- 
ter ellipticals, squares for E/SO galaxies in the field or 
in groups, and triangles for field spiral galaxies. The 
numbers of blue GCs for the survey galaxies were com- 
puted as outlined in Section 12.11 The uncertainties on 
our Tbiue values include Poisson errors on the number 
of GCs and contaminating objects, errors related to the 
magnitude completeness of the data, and errors on the 
galaxy magnitudes and distances (see RZ04 for example 
error calculations). We computed Tbiue for galaxies in 
other studies by combining the published Nqc and blue 
fraction with the galaxy absolute magnitude and appro- 
priate M/Lv- The error bars are based on the published 
errors on Nqc along with uncertainties on the galaxy 
magnitude and distance. Quantities shown in Figure ^ 
are listed in Tabled 

The data in Figure ^ show an overall trend of increas- 
ing Tbiuo with increasing galaxy mass. The two most 
massive galaxies included — the Virgo cluster ellipticals 
NGC 4406 and NGC 4472 — have 2-2.5 times the Tbiue 
values of the lower-mass field galaxies on the left side 
of the plot. A factor that may contribute to this trend 
bears discussion before wc attempt to interpret the data. 
Following ZcDf & Ashman (199^, we used a constant 
M/Ly to calculate Mq for the elliptical galaxies. In 
fact, the stellar mass-to- light ratios of ellipticals almost 
certainly vary with galaxy luminosity, because more lu- 
minous ellipticals are redder (this is the well-known color- 
magnitude relation). The redder stellar populations of 
more luminous ellipticals invariably have higher M/Ly 
in any standard stellar populations model. 

The detailed dependence of M/Ly with luminosity 
goes as M/Ly oc L ^^-^"^ if the color-magnitude relation 
is due to metallicity effects (e.g.. Dressier et al. 1987), 
or as M/Ly oc L if most of the effect is due to 

age differences (Zepf & Silk 1996 and references therein). 
If we assume the latter variation of L°-^°, then Tbiuc 



TABLE 1 
Properties of Galaxies in Figure 



Name 


Type 




Mass 


Tbiuc 


Reference 






[log(M/M0)] 






NGG4472 


E2 


—23.1 


12.2 


2.6±0.4 


5 


NGC4406 


E3 


-22.3 


11.9 


2.5±0.4 


8 


NGC4374 


El 


-22.1 


11.8 


1.9±0.4 


4 


NGC4594 


SO 


-22.4 


11.8 


2.0±0.3 


8 


NGC5128 


EG 


-21.9 


11.6 


1.1±0.2 


5 


NGC3379 


El 


-20.9 


11.3 


1.0±0.2 


8 


NGC1052 


E4 


-21.0 


11.3 


0.9±0.3 


3 


NGC7814 


Sab 


-20.4 


10.9 


1.3±0.5 


7 


NGC4157 


Sb 


-20.4 


10.9 


0.6±0.2 


9 


NGC2683 


Sb 


-20.5 


10.9 


0.9±0.3 


9 


M31 


Sb 


-21.8 


11.4 


1.1±0.2 


1, 2 


Galaxy 


Sbc 


-21.3 


11.2 


0.9±0.2 


1, 10 


NGC3556 


Sc 


-21.2 


11.1 


1.4±0.5 


9 



References. — (1) Ashman & Zepf 1998; (2) Barmby et 
al. 2000; (3) Forbes, Georgakakis, & Brodie 2001; (4) Gomez & 
Richtler 2004; (5) Harris, Harris, & Geisler 2004; (6) Rhode & 
Zepf 2001; (7) Rhode & Zepf 2003; (8) Rhode & Zepf 2004; (9) 
Rhode & Zepf 2005; (10) Zinn 1985. 



Note. — NGC4472, NGC4406 and NGC4374 are Virgo cluster 
galaxies; all others are in the field or in groups. My was derived as 
follows: for NGC4472, by combining V from RZOl with m—M from 
Whitmore et al. 1995; for NGC4406 and NGC4374, combining V 
from RC3 (deVaucouleurs et al. 1991) with m — M from Whitmore 
et al. 1995; for NGC3556, from RC3; for NGC5128, combining V 
from Dufour et al. 1979 with m — M from Israel 1998; for the 
Milky Way and M31, from Ashman & Zepf 1998; and for all other 
galaxies, combining V from RC3 with m — M from Tonry et al. 
2001. 



for the luminous Virgo ellipticals would be reduced by 
a factor of ~1.2 relative to the value for the less mas- 
sive elliptical, NGC 3379. Since Tbiuc for the massive 
Es is more than twice that of the lower-mass galaxies in 
the figure, possible M/Ly variations can only account 
for about one-third of the observed difference in Tbiuo- 
Thus the trend of increasing Tbiuo with increasing host 
galaxy mass is significant and is likely to have important 
implications for galaxy and GC formation. 

2.4. Interpretation and Implications for Galaxy 
Formation 

The prediction of AZ92 that spiral and elliptical galax- 
ies should have similar relative numbers of metal-poor 
GCs can be directly tested using the data in Figure ^ 
Tbiuo for the six spiral galaxies has a weighted mean value 
of 1.2±0.1. Tbiue for the massive Virgo cluster ellipticals 
is twice that of the spirals. We noted earlier that, even 
considering the possible effects of stellar populations vari- 
ations among ellipticals, the difference in Tbiuo between 
the spirals and cluster ellipticals is significant and indi- 
cates that, in terms of the total numbers, combining the 
metal-poor GC populations in a typical spiral galaxy is 
not sufficient to produce the metal-poor GC population 
of a luminous elliptical galaxy. A different — or addi- 
tional — formation mechanism besides the straightfor- 
ward merging of AZ92 must be responsible for creating 
the populous metal-poor GC populations of luminous el- 
lipticals in high-density environments like those plotted 
with filled circles in Figure 1. 

More generally, the data in Figure 1 show a rough trend 
of larger relative numbers of metal-poor GCs with in- 
creasing galaxy mass. Such a trend could be explained 
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by a scenario in which galaxies are formed hierarchically 
and metal-poor GC formation occurs over a finite period 
in the early Universe. As we noted earlier, in hierarchical 
pictures more massive galaxies have a greater fraction of 
their mass assembled by a given redshift than less mas- 
sive galaxies. If GC formation was favored at a specific 
epoch early in the Universe's history, then more massive 
galaxies would consequently have a greater fraction of 
their mass in GCs. Because the GCs that formed in this 
early epoch would be metal-poor, the result would be 
that massive ellipticals hke NGC 4472 and NGC 4406 in 
the Virgo cluster would today have higher Tbiuc values 
than less massive galaxies. In this type of picture, lower- 
mass ellipticals like NGC 3379 start their formation later 
and therefore form a greater fraction of their stars at a 
later time, when metal-poor GC formation is no longer 
favored. Thus they have lower mass-normalized specific 
frequencies of metal-poor GCs. 

Santos (2003) suggested that the type of scenario just 
outlined is indeed responsible for differences in the blue 
GC specific frequencies of galaxies of different masses 
and environments. In the proposed picture, metal-poor 
GCs formed at high redshift when gas-rich protogalac- 
tic building blocks merged to create larger structures. 
This structure formation is suppressed temporarily by 
reionization. Stellar evolution continues to enrich the in- 
tergalactic medium in the intervening time, so that once 
structure formation resumes, any GCs that are created 
are comparatively metal-rich. Because in this picture 
metal-poor GC formation occurs over a finite period, and 
the redshift at which a given galaxy began its assembly 
varies according to the final mass of the galaxy and its 
location in the Universe, differences in Tbiuo are produced 
for galaxies of different masses and environments. 

Detailed modeling of how the relative numbers of blue 
GCs should change with galaxy mass and environment 
for a given GC formation epoch has yet to be accom- 
plished. We have, however , compared ou r data with a 
preliminary calculation from lSantos 1 1|2003() in which blue 
GCs form before z = 10. The Tbiuc differences we ob- 
serve are no larger than what would be expected if metal- 
poor GCs are formed prior to this redshift. Furthermore, 
the observed differences appear to be smaller than what 
is expected if metal-poor GCs form at redshifts much 
greater than 10. More theoretical work is required be- 
fore detailed comparisons between data and model pre- 
dictions can be made. But the idea of a preferred early 
formation epoch for the metal-poor GCs, combined with 
biased hierarchical galaxy formation, is appealing as a 
possible explanation for the apparent trend in the data. 

As described in Santos (2003) , the observed spatial dis- 
tributions of metal-poor GCs can help test the idea of a 
biased, hierarchical origin for galaxies and their GC sys- 
tems. This is because structure that forms before a fixed 
redshift threshold is more centrally-concentrated than 
the dark matter distribution of the final galaxy. Like- 
wise GCs formed within higher-mass protogalactic frag- 
ments will be more concentrated toward the center of a 
galaxy's dark matter halo than those formed in lower- 
mass fragments. Testing whether specific hierarchical 



models work requires predictions for the metal-poor GC 
spatial distributions of galaxies with a range of masses 
and environments for comparison with our observed ra- 
dial distributions (RZOl, RZ03, RZ04). 

Perhaps relevant to this discussion is a recent study 
by Strader, Brodie, & Forbes (2004) that finds a weak 
correlation between the location of the blue GC peak in 
a sample of giant galaxies and the galaxies' luminosities. 
They state that the presence of the correlation means 
that the metal-poor GCs "knew" about the galaxy to 
which they would eventually belong and may have been 
located within that galaxy's dark matter halo when they 
formed. We suggest that within the hierarchical frame- 
work outlined above, one might expect that because the 
assembly process begins earlier for luminous galaxies in 
denser environments, chemical enrichment (for similar 
reasons) proceeds more quickly and efficiently in these 
locations, so that metal-poor GCs formed in these galax- 
ies are slightly more metal-rich than their counterparts 
in less luminous field galaxies. This could produce a shal- 
low positive correlation between the mean metallicity of 
the metal-poor GCs and galaxy total luminosity. Fur- 
ther study of how the blue GC systems of galaxies in the 
field compare to those of galaxies in dense environments 
would help to explore this possible explanation. 

We have emphasized here the links that the metal-poor 
GCs have with galaxy formation, but the metal-rich GCs 
are perhaps equally important as signposts of what has 
occurred over the course of a galaxy's history. In the 
type of hierarchical picture described, the galaxy assem- 
bly process begun at high redshift would be expected to 
continue to the current epoch, with the giant elliptical 
and spiral galaxies of today being the product of many 
major and minor mergers. Formation of metal-rich GCs 
was predicted and has been observed to occur during 
these mergers, provided they are gas-rich (e.g., AZ92, 
Schweizer 1987, and many others). The proportion of red 
and blue GCs, as well as how broad the metallicity range 
of the final metal-rich GC peak is, would then depend on 
the exact merger and assembly history of a particular 
galaxy. This idea is borne out in a general sense by re- 
sults from semi-analytic galaxy formation models (e.g., 
Beasley et al. 2002), which give rise to broadband GC 
color distributions that have a large range in appearance 
and that refiect the wide variation of possible galaxy ac- 
cretion histories. The GC color distributions we observe 
likewise exhibit a range of morphologies (RZOl, RZG4). 
Our survey data may also suggest that Tied, the number 
of red GCs normalized by galaxy stellar mass, increases 
with galaxy mass (RZ04). This might be accounted for 
in hierarchical merger models if a larger number of gas- 
rich mergers contribute to the assembly of more massive 
ellipticals located in richer environments. 
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